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Abstract We study the establishment of three-planet resonances - similar to the 
Laplace resonance in the Galilean satellites — and their effects on the mutual in- 
clinations of the orbital planes of the planets, assuming that the latter undergo 
migration in a gaseous disc. In particular, we examine the resonance relations that 
occur, by varying the physical and initial orbital parameters of the planets (mass, 
initial semi-major axis and eccentricity) as well as the parameters of the migra- 
tion forces (migration rate and eccentricity damping rate) , which are modeled here 

M— I , through a simplified analytic prescription. We find that, in general, for planetary 

masses below 1.5 Mj, multiple-planet resonances of the form n^:n-2\n\=\:'l-A and 
1:3:6 are established, as the inner planets, m\ and m2, get trapped in a 1:2 reso- 
J^ nance and the outer planet 7713 subsequently is captured in a 1:2 or 1:3 resonance 

with mi. For mild eccentricity damping, the resonance pumps the eccentricities 
of all planets on a relatively short time-scale, to the point where they enter an 
inclination-type resonance (as in Libert & Tsiganis 2011); then mutual inclina- 
tions can grow to ~ 35°, thus forming a "3-D system". On the other hand, we find 
that trapping of 7712 in a 2:3 resonance with mi occurs very rarely, for the range 
^. \ of masses used here, so only two cases of capture in a respective three-planet res- 

onance were found. Our results suggest that trapping in a three-planet resonance 
can be common in exoplanetary systems, provided that the planets are not very 
massive. Inclination pumping could then occur relatively fast, provided that ec- 
centricity damping is not very efficient so that at least one of the inner planets 

CO , acquires an orbital eccentricity higher than e = 0.3. 
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1 Introduction 

Recently, the possibility that extrasolar planetary systems can be '3-D systems', 
namely they are composed of two or more planets whose orbital planes have sub- 
stantial values of mutual inclinations, has been considered. Some analytical stud- 
ies have highlighted that such systems can be long-term stable, either following 
normal secular dynamics or due to the action of a phase-protection mechanism, 
such as a mean-motion resonance (MMR) or a secular, Kozai-type resonance (e.g. 
Michtchenko et al. 2006b, Libert & Henrard 2007, 2008, Libert & Tsiganis 2009a). 
A first observational confirmation (McArthur et al. 2010) for the v Andromedae 
system estimated the mutual inclination of the orbital planes of planets c and d of 
this system to ~ 30°. 

In the meantime, a number of studies on the possible formation mechanisms 
of 3-D planetary systems have been realized. Strong dynamical instability of sys- 
tems with giant planets (planet-planet scattering) is generally invoked to explain 
the orbital properties of extrasolar systems, in particular the highly eccentric or- 
bits seen in many (see e.g. Ford, Lystad & Rasio 2005, Fabrycky & Tremaine 
2007, Ford & Rasio 2008, Nagasawa, Ida & Bessho 2008). Many of these have 
shown that the most common outcome of gravitational scattering by close en- 
counters is the hyperbolic ejection of one planet, the 'survivors' having significant 
values of eccentricity and mutual inclination (e.g. Marzari & Weidenschilling 2002, 
Chatterjee et al. 2008, J uric & Tremaine 2008). These works focused on gas- free 
systems, assuming that systems that are relatively compact when the gas nebula 
dissipates will undergo planet-planet scattering after a relatively short time. Only 
a few works investigated the combined action of disc torques and planet-planet 
scattering (e.g. Adams & Laughlin 2003, Moorhead & Adams 2005, Thommes, 
Matsumura & Rasio 2008, Matsumura et al. 2010). However, an unstable crowded 
system - the starting point of typical planet-planet scattering simulations - is not 
the unique result of formation in a gas-disc, since resonant interactions during the 
gas phase should lead to stable resonant systems, such as the ones observed (e.g. 
Michtchenko et al. 2006a, Hadjidemetriou 2006, 2009 for studies of mean-motion 
resonant systems). 

Resonant two-planet systems can also be '3-D', as first shown by Thommes & 
Lissauer (2003) , who studied the evolution of planets trapped in a 1:2 MMR, under 
the effects of gas-driven (Type II) migration, inside the protoplanetary gas disc. 
Libert & Tsiganis (2009b) extended the study of this 'inclination-type resonance' 
to capture in higher order resonances (such as the 2:5, 1:3, 1:4 and 1:5), showing 
that inclination-excitation occurs, as long as the eccentricity of one planet exceeds 
~ 0.4 and the inner planet is not very heavy. 

Concerning resonant systems with more than two planets, it has been recently 
shown that, at least in the case of low-mass planets such as the ones in our solar 
system, gas-driven migration can force the planets to enter into a multiple-planet 
resonance, an analogue of the Laplace resonance in the Galilean moons (Mor- 
bidelli et al. 2007). At present, at least three extrasolar systems are believed to be 
in multiple-planet resonance: HR8799 (e.g. Reidemeister et al. 2009) and G1876 
(Rivera et al. 2010) in a Laplace configuration, and the KOI500 five-planet ex- 
osystem that shows two three-body resonances (Lissauer et al. 2011). Moreover, 
simulations on a possible trapping of HD82943 exosystem in a Laplace resonance 
were also performed by Beauge et al. (2008). Libert & Tsiganis (2011) showed 
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that a Laplace three-planet resonance can be achieved for Jupiter-mass planets, 
but generally becomes unstable as the resonance excites the eccentricities of all 
planets and planet-planet scattering sets in. This mechanism, combining the ac- 
tion of disc torques and planet-planet scattering, typically results in the ejection 
of one of the three planets, leaving behind a two-planet system in a nearly hier- 
archical configuration and with median mutual inclination of ~ 30°. In 10% of 
the simulations presented in that paper, the triple resonance remains stable for 
sufficiently long times, leading to the formation of a '3-D' system, through reso- 
nant inclination excitation, similarly to the two-body inclination-type resonance 
mentioned before. 

As also noted in Morbidelli et al. (2007), a multiple-planet resonance is a 
delicate dynamical configuration; not all resonant ratios can be reached by all 
planetary masses (or mass ratios) and not all resonances are long-term stable. In 
Libert & Tsiganis (2011), we only investigated the dynamical disruption of the 
Laplace 1:2:4 configuration. The goal of the present contribution is to examine 
for what masses of the planets and parameters of the disc a triple resonance can 
be established, and of what type. Also, we examine the frequency of inclination 
excitation, in three-planet resonances. As expected the initial orbital separation of 
the three planets plays a major role in the final outcome and this issue is adressed 
also in the following. 

The paper is organized as follows. In Section 2, we describe the set-up of 
our numerical experiments, while the results of these simulations are presented in 
detail in Section 3. The eccentricity distribution leading to inclination excitation 
is presented in Section 4. Finally, our conclusions are given in Section 5. 



2 Numerical model 

A commonly accepted scenario for the origin of resonant configurations assumes 
disc-induced differential orbital migration of the planets, initially much further 
apart, towards their parent star. While a Jupiter-sized planet carves a gap in the 
gaseous disc, repelling material away from its neighborhood, the massive gaseous 
disc tends to refill the gap, by viscous diffusion. The balance between these pro- 
cesses leads to a relatively slow inward migration of the planets, called Type II 
migration. An order-of-magnitude estimate for the migration rate of a planet was 
given by Ward (1997): 

-9.4«(f)V> (1) 

where P is the orbital period of the body. In the following, the Shakura-Sunyaev 
viscosity parameter is set to a = 4 x 10 -3 (Shakura & Sunyaev 1973) and the 
aspect ratio of the disc to H/a = 0.05. Equation (1) is used to derive an estimate 
for the constant drift rate that we impose on the planet, as in Libert & Tsiganis 
(2009b): we assume a in Equation ([1]) to be constant and equal to the initial value 
of semi-major axis of the planet's orbit. This is of course a simplified prescription, 
but adequate for analyzing the dynamics of resonance encounters. As planet-disc 
interactions also affect the eccentricity of the migrating planet (see Goldreich & 
Tremaine 1980, Papaloizou, Nelson & Masset 2001), exponential damping is also 
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assumed in the following way: 

e 



= -K 
e 



(2) 



where K is chosen to be 0, 1 and 5 in the following. 

We start our simulations with a system of three planets (mi, 7712 and 7713; mx 
being the inner body and 7713 the outer one) evolving around a 1 Mq star. Assum- 
ing that the inner disc is likely to be largely depleted by the time our planets were 
formed, we simulate planet-disc interaction by applying a suitable Stokes-type drag 
force (following Beauge, Michtchenko & Ferraz-Mello 2005) on rri2 and 7713 only. To 
integrate this four-body problem we use a variant of the SYMBA code (Duncan 
et al. 1998), which is able to deal symplectically with close encounters between 
massive bodies, in which the drag terms are added in the equations of motion, 
introducing an exponential drift in semi-major axis (migration) and eccentricity 
(damping). Let us note that, even if the drift is not applied on the inner body, 
when capture in resonance occurs, the two bodies subsequently migrate as a pair 
towards the star - at a slower rate - and both their eccentricities are damped. How- 
ever, for mild eccentricity damping, the growth of the eccentricities may be very 
rapid, as migration in resonance continues (e.g. Lee & Peale 2002, Ferraz-Mello 
et al. 2003). Note that, as shown in Morbidelli & Crida (2007), in more realistic 
simulations of gas-planet interactions the resonance pair may stop migrating (or 
even start migrating outwards), depending on the planetary masses, a and H/a. 
However, in such a case, 7713 would approach a resonance with 777,2 more easily and 
a triple resonance could again be established. To take into account (partly) these 
uncertainties that are related to our simplified migration model, we perform our 
simulations also using migration rates 5 times larger and 5 times smaller than the 
values given by Equation (1). 

In Libert & Tsiganis (2011), the establishment of the Laplace triple resonance 
was reached by a two-step formation process, following the lines of Morbidelli et 
al. (2007): first we simulated the 1:2 resonant capture of the two inner planets, m\ 
and m2 . Then, planet 7713 was introduced in the simulation and forced to migrate 
into the 1:2 MMR with 7712. The Laplace resonance was thus approached, but not 
always established; the orbits of heavy giant planets usually become chaotic and 
planet-planet scattering dissolves the system before the establishment of the 1:2:4 
resonance. Here, the methodology is different, as we do not focus on this particular 
resonance. Instead, we try to identify which triple resonances can be established 
by disc-induced differential orbital migration of three planets. 

Two initial configurations of the three-planet systems are examined in this 
work, depending on the position of the intermediate planet (7712) with respect to 
its 1:2 MMR with m\ (initially located at 1). The semi-major axis of 7712 is set 
to 1.9 (exterior to their 1:2 MMR - see Section [3. ip and 1.4 (interior to their 1:2 
MMR - see Section l3.3[) . In either case, we consider several values of the semi- 
major axis of 7773: 2.8, 3.1 and 3.6 for the exterior case, 2, 2.4 and 2.7 for the 
interior case. The initial eccentricities and inclinations of all planets are nearly 
zero (e = 0.001, i = 0.01°). The effects of initially eccentric orbits are studied in 
Section 13.21 Nine mass configurations are considered for the planets: mi is fixed 
to 1.5 Mj (Mj = Jupiter's mass), while the mass of the second planet is set to 
7772 = 0.75, 1.5 or 3 Mj up and 7773 to 7712/2, 7712 or 2rri2. Our system of units is 
such that G = 1 and Mq = 1. Taking the semi-major axis of mi as the unit of 
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Table 1 MMR captures for three-planet systems of different masses (mi = 1.5 Mj), for 
different values of the eccentricity damping rate. The symbol '*' indicates that the inclination- 
type resonance occurs as the planets migrate in resonance. The migration rate is equal to 
the value given by Equation |[TJ, namely |a'2/a2| = 5.71 X 10 — 6 time units -1 and |a'3/a3| = 
3.19 X 1CT 6 , 2.74 X 10~ 6 or 2.19 X 1CT 6 time units" 1 for a 3 = 2.8, 3.1 and 3.6 respectively. 

a 3 = 2.8 3T 3fi 

m 2 m 3 K=0 K=l K=5 K=0 K=l K=5 K=0 K=l K=5 
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distance (ai = 1), the time unit is defined such that the period of mi is Pi = 2ty. 
All simulations spanned a time equal to 10 6 time units, a time long enough to be 
comparable to the disc's lifetime. The results of these simulations are described in 
the following section. 



3 Results 

3.1 Exterior configuration 

We consider first the evolution in the gas-disc of a three-planet system in which 
the intermediate planet (1712) is initially set exterior to the location of its 1:2 MMR 
with mi: here we use a\ = 1 and a 2 = 1-9. Three different initial positions of the 
outer planet (7713) are used in the following simulations: 2.8 (initially outside the 
2:3 MMR with m 2 ), 3.1 (initially outside the 1:2 MMR with m 2 ) and 3.6 (initially 
outside the 2:5 MMR with m 2 ). For the values of the mass ratios and eccentricity 
damping discussed in the previous section, we aim to determine the frequency of 
(i) capture in different multiple-planet resonance - if capture occurred - and (ii) 
subsequent resonant inclination excitation. 

The results of these simulations are presented in Table [TJ For each mass ratio 
and eccentricity damping value, we describe the time evolution of the system by 
denoting the MMR in which it is captured. Three different outcomes are observed 
in the simulations: 

(a) the unstable case, represented by an empty slot in the Table. The system is 
destabilized in less than 10 time units, without having previously reached any res- 
onant configuratiorlj. As seen in the Table, this occurs primarily for mass ratios 
m 2 /mi > 1 and 7713/mi > 1, i.e. for more massive outer planets, and small orbital 
separations. As K increases, the instability seems to occur only for m 2 /mi > 2 
and 7713/7711 > 2. This implies that, for a massive system in a compact initial con- 
figuration, the eccentricities grow fast and the planetary orbits cross each other, 
before a resonance can be established. 



1 We note that the simulation is stopped if one of the planets is ejected from the system or 
two planets merge. 
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(b) the two-body resonant case. The two inner planets are trapped in a 1:2 res- 
onant configuration, which is found to be destabilized when the outer planet is 
approaching a resonance with 7712. Note that this occurs for an initial a^ = 3.1 or 
higher. For an initial 03 = 2.8, no similar capture of the two inner planets in a 2:3 
MMR is observed. 

(c) the three-body resonant case, where all three planets are trapped in a resonant 
configuration. According to the Table, there are two main resonances, depending 
on the initial configuration (i.e. period ratio P3/P2). Inclination excitation is ob- 
served for both resonances, but only for K < 1. This implies that the excitation 
is related to eccentricity growth in the resonance, as was also shown in Libert & 
Tsiganis (2011). 

As Table [T] shows, massive planets (first lines of Table[T) and an initial location 
of the intermediate planet very close to its 1:2 MMR with the outer planet (03 = 
3.1 - columns 6, 7 and 8 of Table [l| prevent the formation of a multi-resonant 
state. For massive planets, the typical outcome is the ejection of one of the bodies 
from the system at the beginning of the simulation, especially for the particularly 
compact initial configuration a\ = 1, a^ = 1.9 and 03 = 2.8. For the two other 
initial configurations, i.e. larger distance of 7713 from the inner planets, massive 
planets meet easily the 1:2 MMR between mi and 7712, but no capture occurs for 
7713. These two configurations excepted (massive planets and 03 = 3.1), most of 
the simulations result in three-planet resonant configurations. The two dominant 
relations found are: the Laplace 1:2:4 resonance (for an initial an, = 2.8) and 
the 1:3:6 resonance (for an initial 03 > 3). The notation used here means that 
m\ and 7712 are in a 1:2 MMR (712:711 = 1:2), while 7712 and 7713 in a 1:3 MMR 
(713:712=1:3). Thus, the multiple-planet resonance is labeled as 713:772:711=1:3:6. 
The same convention will always be adopted in the following. 

Note that these results may depend on the assumed simplified migration model, 
as well as on the assumed migration rates (see Section f3 . 2 1> . However, we believe 
that the observed trend, i.e. easier establishment of a triple resonance for less 
massive triplets should hold, unless the relative migration of planets is strongly 
affected by variations in the surface density profile of the gas disc (e.g. as in 
Morbidelli & Crida 2007). 

As the system remains in MMR and continues to migrate, the eccentricities 
of the planets become high enough that the system enters an inclination-type 
resonance, which induces rapid growth of the inclinations. Inclination excitation is 
indicated by the presence of the symbol '*' in Table [1] Inclination excitation can 
arise from two resonant mechanisms: the one described by Thommes & Lissauer 
(2003) for two migrating planets in the 1:2 MMR and a similar one acting with 
three planets in a triply-resonant configuration, described in Libert & Tsiganis 
(2011). Examples of both outcomes are shown in Figures [1] and [2] respectively. 

Figure [1] shows the evolution of a three-planet system where two planets go 
through the Thommes & Lissauer's mechanism. While the outer planets migrate, 
mi and 77i2 are captured in a 1:2 MMR at about 3 x 10 time units, characterized by 
the librationof the resonant angles 8\ = Ai— 2A2+TO1 and #2 = X±— 2A2+1372 around 
0°. As a result, the difference of the longitudes of pericenter Aw = vji~vj2 = B1—B2 
starts to oscillate around 0°, i.e. the planets are in apsidal alignment. As both 
planets continue to migrate while in resonance, their eccentricities increase. When 
their values are high enough, the system enters an inclination-type resonance: 
the angles 8^ = 2Ai — 4A2 + 2i?i and Op = 2Ai — 4A2 + 2.O2 start to librate 
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Fig. 1 Inclination excitation for a 1:2 MMR between mi and r»2. The planetary masses are 
ITU = 3 and mo, = 1.5Mj. Initial location of the outer planet is a-j = 3.6. No eccentricity 
damping is considered (K = 0). 



at about 2 x 10 5 time units. A rapid growth of the inclinations of both planets 
is then observed, as well as a libration of the relative longitude of the nodes, 
Afl = Q\ — Q-2. = (6^2 — 9^) /2, around the anti-alignment state. Note that the 
eccentricity and inclination of the outer planet 7713, which is not in resonance 
with the other ones, remain almost constant (i.e. nearly zero) during the whole 
simulation of Figure[T] The body 7773 goes successively through 1:3 (when a^ ~ 1.6 
and a-;, — 3.3) and 1:4 (when a-2 — 0.93 and a-j, — 2.3) commensurabilities with 777,2, 
and at ~ 3.5 x 10 5 time units (not shown on the figure), while approaching 1:5 
commensurability (02 ~ 0.6 and 0,3 ~ 1-76), destabilizes the whole system, leading 
to planet-planet scattering and the subsequent formation of a non-coplanar two- 
planet system, the less massive planet being ejected from the system. Indeed, 
the given hi and masses result in a given change in energy of the resonant pair 
which could be translated to a mean migration rate. If this rate is smaller than 
<X3, then we have divergent migration of 7773. A further analysis of the reason for 
instability in all studied systems (i.e. MMR crossing during divergent migration 
or instability before the resonance) is reserved for future work. Let us note that 
for the same initial configuration of the system of Figure [1] but for less massive 
planets, a capture in a multiple-resonance generally happens, as shown in Table [1] 
(lines 4 to 9). 

Inclination excitation due to multiple-planet resonance is described in Figure 
[2] While the outer planet migrates, the two inner planets are captured in a 1:2 
MMR at about 3 x 10 4 time units, characterized by the libration of both resonant 
angles 6± = Ai — 2A2 + w\ and 82 = \i — 2A2 + zu2 around 0°. From then on, 
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Fig. 2 Inclination excitation for a three-planet system in a Laplace 1:2:4 resonant configura- 
tion. The planetary masses are m.2 = 1.5, m^ = 0.75Mj. Initial location of the outer planet is 
113 = 2.8 and eccentricity damping is set to K = 1. 



these planets migrate as a pair, while the outer planet is captured in a 1:2 MMR 
with the second body at about 7 x 10 4 time units (one of the resonant angles 63 = 
\2-2\3-\--cd2 and 84 = A 2 — 2A3+ci73 is in libration). Thus the system is captured in 
a Laplace-type resonance, whose critical angle is <f> = \\ — 3A 2 + 2A3Q As the three 
planets continue to migrate while in resonance, their eccentricities increase. When 
their values are high enough, the system enters an inclination-type resonance: the 
angles 6p = 2Ai - 4A 2 + 2Qi, Op = 2Ai - 4A 2 + 20 2 , 6^2 = 2A 2 - 4A 3 + 2Q 2 and 
9p = 2A 2 — 4A3 + 2L?3 start to librate around 180° at about 3 x 10 J time units, 
and a rapid growth of the inclinations of the three planets is observed. 

The same type of evolution is observed for higher-order multiple-planet reso- 
nances. Figure [3] illustrates the inclination excitation for a three-planet system in 
a 1:3:6 resonant configuration (mi and m 2 are in a 1:2 MMR, m 2 and 7713 in a 1:3 
MMR), whose resonant angle is <j> = 2Ai — 5A 2 + 3A3. The initial configuration of 
the system is the same as the one of Figure Q] but the planets m 2 and iris are less 
massive, which favors a 1:3 resonance capture of TO3 during its migration towards 
m 2 . One of our simulations even ends in a 1:4:8 resonance capture (mi and m 2 
are in a 1:2 MMR, m 2 and m^ in a 1:4 MMR - see Figure[3]), characterized by the 
libration of <f> = 3Ai — 7A 2 +4A3. However, the inclination excitation mechanism is 
not effective in this case, as the system is destabilized shortly after the beginning 
of the libration of the inclination-type resonant angles. 



2 The critical angle is a combination of the resonant angles 8j, which must respect the 
d'Alembert characteristics; thus, the sum of all integer coefficients that appear in front of each 
mean longitude should be zero. 
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Fig. 3 Inclination excitation for a three-planet system in a 1:3:6 resonant configuration. The 
planetary masses are m,2 = 1.5, m.3 = O.ToMj. Initial location of the outer planet is 0,3 = 3.6. 
No eccentricity damping is considered (K=0). 
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Fig. 4 Inclination excitation for a three-planet system in a 1:4:8 resonant configuration. The 
planetary masses are ni2 = 3, m-j, = 1.5Mj. Initial location of the outer planet is 03 = 3.6 and 
eccentricity damping is set to K = 1. 
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Fig. 5 Left: Effect of eccentricity damping of the form e/e = —K\a/a\ with K = 1. The 
planetary masses are 7712 = 0.75, m^ = 0.75Mj. Initial location of the outer planet is 0,3 = 2.8. 
Inclination excitation takes place at ~ 1.6 X 10° time units, before the disruption of the system 
at ~ 2.1 X 10° time units. Right: Same as left Figure, but for K = 5. Due to the high value 
of the eccentricity damping, the growth in eccentricity is not sufficient to produce inclination 
excitation, so that inclinations remain bounded to tiny values (less than 0.02°). 



Resonant inclination excitation appears to be quite common in our model, as 
long as eccentricity damping is not too strong. Indeed, a damping rate of K = 5 (see 
Equation ([2])) does not allow eccentricities to increase sufficiently for an inclination- 
type resonance to occur, as shown by the example of Figure[5] A similar feature has 
also been highlighted in the case of a two-planet system in Thommes & Lissauer 
(2003) and Libert & Tsiganis (2009b). Lee & Thommes (2009) studied the capture 
of a two-planet system in the 1:2 MMR and found that inclination excitation occurs 
only for K ~ 1 or less, for a range of mass ratios and migration rates similar to 
what we adopt here. Our simulations confirm this limit and show that it applies 
also to the case of three-planet resonancqj. 

For no or mild eccentricity damping (i.e. K < 1), Table [1] suggests that 65% 
of the triply-resonant systems can go through inclination-type resonance. Hence, 
significant mutual inclinations between the orbital planes of a three-planet system 
seems to be a common feature, in the formation scenario depicted by our model. 

Mechanisms responsible for halting Type II migration are not well understood 
so far. In the case of mild (or no) eccentricity damping, the triply resonant con- 
figuration typically becomes unstable, since the eccentricities of the planets keep 
increasing as migration continues. However, if the disc dissipates before the system 



3 In that paper another, faster, mode of inclination-excitation was found, not clearly asso- 
ciated with the libration of any critical angle. We have not found an analogue to this mode 
within our results 
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Fig. 6 Evolution of the system shown in Figure \3\ after suddenly switching off the migration 
force (at 4 X 10 5 time units). The system remains stable, with high mutual inclinations values 
between the orbital planes. 



dissolves, a stable configuration can be reached. If it occurs during the inclination- 
type resonance phase, it can produce a stable resonant non-coplanar three-planet 
system. An example of such a system is given in Figure [6l where the migration of 
the outer planets is abruptly switched off at 4 x 10 time units in the simulation of 
Figure [3] (1:3:6 resonant system). The system remains stable, in orbital resonance 
(see the constant evolution of the semi-major axes), with small variation in its 
eccentricities, and with mutual inclinations of the orbital planes between 5° and 
35°. 

An obvious limitation of our model is that we do not account for possible in- 
clination damping by the gaseous disc. Inclined planets would spend most of their 
orbit away from the disc plane, so it is not easy to find empirical formulae for incli- 
nation damping by performing hydrodynamical simulations. On the other hand, if 
eccentricity damping in discs is moderate (a requirement for the inclination-type 
resonances to be reached), one should expect the same for inclination damping. 



3.2 Varying the initial eccentricities and migration rates 



In the previous section, we assumed planets to be initially on (nearly) circular 
orbits. The effects of initially eccentric orbits were also studied, in order to check 
whether capture into multiple-planet resonance can occur for non-circular orbits. 
The results are given in Table [2j Three different initial states were analyzed: (A) 
d = 0.15, e 2 = 0.001, e 3 = 0.001, (B) e a = 0.001, e 2 = 0.15, e 3 = 0.001 and 
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Table 2 Same as Table ^ for different values of the initial eccentricities of the planets. A, 
B and C denote different initial configurations: (A) e\ = 0.15, e 2 = 0.001, e 3 = 0.001, (B) 
ei = 0.001, e 2 = 0.15, e 3 = 0.001 and (C) ei = 0.001, e 2 = 0.001, e 3 = 0.15. 
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(C) ei = 0.001, e2 = 0.001, e% = 0.15. For each case, we performed the same 
set of simulations as in the previous section. As seen from the Table, the main 
result is that, for no or mild eccentricity damping, the system is easily destabilized 
before the establishment of any MMR, so much fewer multiply-resonant systems 
are formed. Most of the inclination excitation observed is due to the 1:2 MMR 
between the two inner planets. For strong eccentricity damping [K = 5), the 
results are much more similar to those presented in the previous section. 

The new result that can be found in Table [2] is the possibility for planets on 
initially non-circular orbits of being locked in a 1:2:5 resonance (i.e. mi and mi 
are in a 2:5 MMR, m 2 and m 3 in a 1:2 MMR). 

The results presented above are by no means a complete study of the de- 
pendence of resonance capture on the initial planetary eccentricities. This subject 
certainly deserves a more thorough analysis. On the other hand, our results clearly 
show that the probability of capture of three planets in a multiply-resonant, quasi- 
stable state is much lower, when the initial non-resonant orbits are moderately 
eccentric and eccentricity damping is not strong. 

In the previous section, the migration rates of the outer planets were taken 
from Equation ([TJ. As explained in that section, this is a rough approximation 
and therefore the results should be taken with care, as the actual migration mode 
of the planets may be different. Hence, the influence of the migration rate on the 
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Table 3 Same as Tabled for different values of the migration rate of the planets. 
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Fig. 7 Trapping in a 1:2:5 resonance (mi and mi are in a 2:5 MMR, rri2 and ms in a 1:2 
MMR). The planetary masses are m^ = 1.5, 7713 = 1.5Mj. Initial location of the outer planet 
is ag =3.1 and eccentricity damping is set to K = 1. 
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Table 4 Same as TablefTlfor the interior configuration. The migration rate is equal to the value 
given by Equation JTJ, namely |a'2/a2| = 9.03 X 10~ 6 time units -1 and 1 03/03! = 5.29 X 10 — 6 , 
4.02 X 10~ 6 or 3.37 X 10~ 6 time units -1 for a 3 = 2, 2.4 and 2.7 respectively. 
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results found in the previous section has to be examined. We analyzed the behavior 
of the same set of simulations as before, evolving now under five times larger and 
five times smaller migration rates. As can be seen in Table [3l the outcomes of 
the simulations are similar. A more careful examination of it shows that actually 
a slower migration rate favors both the capture in a triple resonance as well as 
inclination excitation; however the difference is not very dramatic, taking into 
account that the range of values of \a/a\ used here corresponds to an order of 
magnitude. Let us note that one simulation differs again from all the others: the 
two inner planets are trapped in a 2:5 MMR, before the establishment of the 
1:2:5 resonant configuration (mi and ni2 are in a 2:5 MMR, m-2 and 1713 in a 1:2 
MMR, see Figure [7])- The resonant angle is 4> = 2Ai — 8A2 + 6A3. The system 
is disrupted before an inclination-type resonance can be established, due to the 
large eccentricities reached; apparently the inclination resonance is located for even 
higher values of e than required for the system to remain stable. The trapping of 
mi and mi in a high order resonance occurs for a slower migration rate; this 
observation was also highlighted in Libert & Tsiganis (2009b). 



3.3 Interior configuration 



In the previous simulations, mi was initially located exterior to its 1:2 MMR 
with m\. This means that, during its orbital migration, mi meets (and is usually 
captured into) this first-order resonance. One may wonder what would happen if 
mi and mi were initially located interior to their 1:2 MMR, or manage to jump 
across the 1:2 MMR during migration, as is thought to have been the case for the 
Jupiter-Saturn system (see Pierens & Nelson 2008). In the following, we set the 
initial semi-major axis of mi to 1.4, and consider three initial values for 03: 2, 2.4 
and 2.7. 

For the initial conditions considered above, we have performed the same set of 
81 simulations as in Table[T] MMR captures are shown in Tabled Less than 10% 
of the simulations result in a two- or three-planet resonance capture. Most of the 
simulations show a quick destabilization of the system (planet-planet scattering or 
even mergers). This is most likely an outcome of the initially compact configuration 
of the system and relatively large assumed masses. In five simulations, mi and mi 
manage to get locked in a 2:3 MMR; this happens when mi < mi. Only two 
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Fig. 8 Inclination excitation for a three-planet system in a 3:8:12 resonant configuration. The 
planetary masses are m,2 = 0.75, mo, = 1.5Mj. Initial location of the outer planet is 03 = 2.7 
and eccentricity damping is set to K = 1. 



simulations, corresponding to mi = 1.5, 7712 = 0.75 and 7713 = 1.5Afj, result in 
three-planet resonant configurations: a 4:10:15 resonance (mi and 7772 are in a 
2:3 MMR, m.2 and 771,3 in a 2:5 MMR), characterized by the libration of <j) = 
6A1 — IIA2 + 5A3, and a 3:8:12 resonance (mi and mo, are in a 2:3 MMR, 771,2 
and 7773 in a 3:8 MMR - see Figure [8]), characterized by the libration of </> = 
10Ai — I8A2 + 8A3. For this last multiple-planet resonance, an inclination-type 
resonance is also observed at large-enough eccentricities (see the libration of the 
resonant angles 9p = 4Ai-6A 2 +2J?i, 9p = 4Ai-6A 2 +2J? 2 , Op = 6A 2 -16A 3 -|-10J?2 
and 0,2 = 6A 2 - I6A3 + 10« 3 ). 

Given the results of this section, it seems that the formation of a three-planet 
resonant configuration where 7711 and 7772 are locked in a 2:3 MMR, is rather 
unlikely, at least for planetary masses larger than ~ 1A1 j and mass ratios 7772/777,1 > 
0.5. On the other hand, it may be common for mass ratios similar to the Jupiter- 
Saturn pair. This topic deserves further investigation in the future. 



4 Eccentricity distribution 



In this section, we discuss the eccentricity distribution leading to inclination ex- 
citation, for the masses, migration rates and initial configurations of the systems 
presented in Tables [1] and [3] (i.e. initially circular orbits). We plot in Figure[9l the 
eccentricity values of the three planets, measured just before the increase of the 
inclinations. All systems of Tables [1] and [3] that lead to inclination excitation are 
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Fig. 9 Eccentricity distribution before inclination-type resonance for all systems of Tables [T] 
and [3] that lead to inclination excitation 



shown in this figure. We see that inclination excitation occurs for a wide range of 
eccentricities that depends strongly on the masses of the two inner planets. For 
mi < 77i2, the eccentricity of mi is higher than the one of 7712, and vice versa. The 
eccentricity of the outer body 7773 is rather small for all mass ratios considered here. 
It is interesting to note that left panel of Figure [9] shows an apparently correlation 
between e^ and ei (see the linear regression). Moreover, no system with eccentric- 
ities e\ and e-2 both smaller than 0.3 can be found on this plot. This means that 
inclination-type resonance requires at least one of the two inner planets to reach 
an eccentricity higher than e = 0.3. 

The eccentricity distribution given in Figure [9] is close to the final eccentricities 
distribution of the resonant systems. Indeed, when the system is in the inclination- 
type resonance, the eccentricities keep growing but slowly (see e.g. Figures 2 and 
3), for as long as the system can remain stable. 



5 Conclusion 



In the present work we studied the formation of three-planet resonances in ex- 
trasolar planetary systems. Such resonances can be established during the phase 
of gas-driven planet migration. If the resonance persists as the gas dissipates, a 
stable resonant three-planet system can form. Otherwise, the system is dynami- 
cally disrupted, typically leading to the ejection of one planet. Libert & Tsiganis 
(2011) studied the 1:2:4 resonance and showed that, if the system becomes unsta- 
ble, the survivors tend to have large orbital separations (i.e. hierarchical systems). 
Of particular interest is the evolution of the orbital inclinations of planets trapped 
in a 1:2:4 resonance; as migration continues and the eccentricities of the planets 
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grow, inclination-type resonances come to play. Thus, similarly to the Thommes 
& Lissauer mechanism (for two planets), the inclinations can grow to values larger 
than ~ 30°. 

Here we showed that trapping to several triple resonances can occur, depending 
of course on the masses of the planets, their initial orbital separations, and the 
parameters of the disc - described here in a simplified way through the values of 
the migration and eccentricity-damping rates. Large planetary masses typically 
lead to instabilities that destroy the system before a resonance can be established. 
For masses smaller than ~ 1.5 Mj, the most common resonant relations found were 
the ri3:n2:ni = l:2:4 and 1:3:6. These resonances appear to be stable for relatively 
long times, especially when eccentricity damping is applied. For mild eccentricity 
damping, these resonances can pump the eccentricities of all planets on relatively 
short time-scales. Then, the planets can enter an inclination-type resonance (as in 
Libert & Tsiganis 2011), which pumps their inclinations up to ~ 35°, thus forming 
a "3-D system". 

When iri2 is started interior to its 1:2 MMR with mi, the strongest resonance 
in which it could be trapped is the 2:3 MMR with mi. This occurred only in very 
few cases and, consequently, only two cases of capture in a respective three-planet 
resonance were found. We note that, according to recent results of hydrodynamical 
simulations (e.g. Pierens & Nelson 2008), trapping in the 2:3 MMR is considered 
to be the most likely result of gas-driven migration for the Jupiter-Saturn pair in 
our solar system. The fact that we do not find such configurations may reflect our 
large planetary masses (and mass ratios), compared to the Jupiter-Saturn system. 
Hence, it may actually be the case that trapping in the 2:3 MMR (and related 
three-planet resonances) is possible only for low-mass planets. 

Our results suggest that trapping in a three-planet resonant configuration dur- 
ing gas-driven migration may be common, as long as the planets in question are 
not very massive. Also, if eccentricity damping is not too strong in discs, resonant 
inclination pumping can be efficient, provided that one of the two inner planets 
eccentricities is higher than 0.3. We note however that our results should be taken 
with care, as our migration model is not very advanced — this is why we have 
avoided presenting statistics on the efficiency of trapping, etc. On the other hand 
we have no reason to question the generic character of the resonance-trapping and 
inclination-excitation mechanisms. In order to be able to evaluate the probability 
of capture in a given triple resonance, its stability time (compared to the disc- 
evaporation time scale) and its efficiency in exciting the inclinations, we need to 
resort to more advanced models of migration - e.g. the one used in Thommes et 
al. (2008) - that take into account variations in the gas surface density induced by 
the relative motion of the planets, as well as gas dissipation. We hope to be able 
to report on these issues in a forthcoming paper. 
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